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In the first part of this contribution we present a numerical study motivated by recent attempts to 
understand the nonperturbative aspects of QCD at temperatures T ~ a few times the deconfine- 
ment temperature T c , We focus on the pure gauge theory, and ask whether the deficit in pressure 
and entropy, with respect to their free-gas values, is particular to SU (3). We find that the deficit 
in SU(4),SU(8) for T < 2T C , l.6T c , respectively, is remarkably close to that of SU(3). This sug- 
gests a similar deficit for SU (°°), which is fortunate since this theory is simpler, and can serve to 
constrain the possible dynamics underlying the deficits. 

In the second part we seek for signs of a Hagedorn temperature 7# in pure lattice SU (N) gauge 
theories with N = 8, 10, 12. Since one expects Th > T c , we measure masses of strings in the 
metastable confined phase above T c , and extrapolate them to zero to estimate Th- For SU (12) we 
find that Th/T c = 1.1 16(9), when we extrapolate with a critical exponent of the three dimensional 
XY model, which seems to be preferred over a mean field exponent by our data. 
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1. Bulk thermodynamics of lattice SU (N) gauge theories at large-TV 

In this section we present a study of bulk thermodynamical properties of 5*7 (4), and 5*7 (8) 
pure lattice gauge theories and compare to a study by Boyd et al. of SU (3) [|T|]. Recent calculations 
of various properties of SU (N) gauge theories [||]-[||] have demonstrated that N = 8 is very close to 
N = oo for most purposes. These also provide information on the critical coupling for various lattice 
sizes and N. Thus our calculations should provide us with an accurate picture of what happens to 
bulk thermodynamics at N = °°. For a more detailed version of this section we refer to [01 . 

1.1 Lattice setup and methodology 

We define the gauge theory on a discretised periodic Euclidean four dimensional space-time 
with L] x L t sites, and perform Monte-Carlo simulations of a simple Wilson action. We use the 
Kennedy-Pendelton heat bath algorithm for the link updates, followed by five over-relaxations of all 
the SU (2) subgroups of SU (N). To evaluate bulk thermodynamics we choose to use the "integral 
method" (see for example [[]]]) 1 and so the pressure p and interaction measure A are 



p/T 4 = 6L 4 I df3'((u p ) T - (u p ) ), A/T 4 = 6L 4 ((u p (P)) - (u p (f3)) T ) x 



9^ 



(1.1) 

Here {u p )j is the plaquette average on a T > lattice with L t < L s , while (u p )o is measured on a 
lattice with relatively large L t = L s . To evaluate the scaling of the temperature T with the coupling 
j8 = 2N/g 2 , and the derivative in Eq. (|j~I|), we use calculations of the string tension, a, in lattice 
units (e.g. [§]). Other physical choices to fix the scale differ in modest 0(a 2 ) differences. Finally 
the value of T c / y/o for the different gauge groups and L ? = 5 in taken from @ 

We performed calculations of (u p )t in SU (4) on 16 3 5 lattices and in SU (8) on 8 3 5 lattices for 
a range of /} values corresponding to T/T c € [0.89, 1.98] for 5*7(4), and to T /T c £ [0.97, 1.57] for 
51/(8). S ince we use Lt — 5, while the data for 5*7 (3) in [jl|] is for Lt — 4, 6, 8, we also performed 
simulations for 5*7(3) on 20 3 5 lattices with T/T c £ [1,2]. The measured plaquette averages are 
presented in 

We performed the 'T = 0' calculations of (u p )o on 20 4 lattices for SU(3), and on 16 4 lattices 
for SU(4), which are known to be effectively at T = [||, ||] for the couplings involved. For SU(8) 
however, using 8 4 lattices would not be adequate for the largest j8 -values, and we take instead the 
SU(8) calculations on larger lattices in [gj, and interpolate between the values of /3 used there with 
the ansatz (w p ) (j3) = {u p )% T (P) + ^^i(aVa) 4 + c 4( g 8 + c 5 g w , where {u p )% T (p) is the lattice 
perturbative result to 0(g 6 ) from [[|] and = 8. Our best fit has # 2 /dof = 0.93 with dof = 2, and 
the best fit parameters are c\ = —6.92, C5 = 26.15, and a gluon condensate of = 0.72. 

1.2 Finite volume effects 

For N = 4, 8, one is able to use lattice volumes much smaller than what one needs for 5*7 (3) 
(like those in |jl]]) as the longest correlation length decreases rapidly with Af [|[ Qj. This is also 
theoretically expected, much more generally, as N — > °°. The main remaining concern has to do 



'Our lattices are too coarse for the "differential" method, and we found that the Wang-Landau algorithm for the 
evaluation of the density of states, did not converge easily for SU (8). 
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with tunneling configurations, which occur only at T c when V — > °°. On our finite volumes, this is 
no longer true, and we minimise finite-V corrections by calculating the average plaquettes only in 
field configurations that are confining, for T < T c , or deconfining, for T > T c . For SU(3), where 
the phase transition is only weakly first order, it is not practical to attempt to separate phases. This 
will smear the apparent variation of the pressure across T c in the case of SU(3). 

To confirm that our finite volume effects are under control we have compared the SU(8) value 
of (u p (P)) as measured in the deconfined phase of the our 8 3 x 5 lattice with other x 5 results 



from other studies [jllj], and find that the results are consistent at the 2 sigma level. We perform a 
similar check for the confined phase on the same lattices, and again find that finite volume effects 
are small, mostly on a one sigma level. The data supporting these checks, together with the checks 
related to the next paragraph, is presented in [^]. 

A similar check for the confined phase on L 4 lattices leads to conclude that a size L = 8 in 
SU(8) is not large enough, as we find that the plaquette average has a significant change (on a 
16 sigma level) from L = 8 to L = 16 for our largest value of /3. By contrast, for St/ (4) the 
finite volume effects seem not to be large on the 16 lattice as we checked for our largest value of 
j8 = 1 1.30. There the value of the plaquette on a 20 4 lattice is consistent within ~ 2.3 sigma with 
the value on a 16 4 lattice. 



1.3 Results 

In presenting our results for the pressure, we shall normalize to the lattice Stephan-Boltzmann 
result given by (/V^ 4 )f ree _ gas = (N 2 — l)fg xRj(Lt)- Here Rj includes the effects of discretization 
errors in the integral method [j^fTj]. The same normalisation is applied for the interaction measure 
A, and for the energy density e = A + 3p, and entropy density s = (A + 4p) /T. 

We present our N = 4 and N = & results for p/T A in the left plot of Fig. [[]. We also show our 
calculations of the SU (3) pressure for L t = 5, as well as the L t = 6 calculations from [jjj]. In the 
right plot of Fig. [T] we present results for the normalized energy density e, and normalized entropy 
density s. The lines are the SU (3) result of [[I]] with L t = 6. Again we see very little dependence on 
the gauge group, implying very similar curves for N = °°. One can clearly infer that the pressure 
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Figure 1: The pressure (on the left), and energy and entropy densities (on the right) normalized to the 
lattice Stephan-Boltzmann pressure, including the full discretization errors given by Rj. In the pressure 
plot, the symbols' vertical sizes are representing the largest error bars (which are received for the highest 
temperature). The lines are the St/ (3) results and L t — 6 from [|l]]. 
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in the SU (4) and SU (8) cases is remarkably close to that in SU (3) and hence that the well-known 
pressure deficit observed in SU(3) is in fact a property of the large-A^ planar theory. This implies 
that the dynamics that drives the deconfined system far from its noninteracting gluon gas limit, 
must remain equally important in the N = °° planar theory. This is encouraging since that limit is 
simpler to approach analytically, for example using gravity duals, and also because it can serve to 
constraint and point to important ingredients of analytical approaches (see the discussion in [0]). 

2. In search of a Hagedorn transition in SU (N) lattice gauge theories at large-Af 

In this section we present a calculation of the change in the mass m(T) of Polyakov loops in 
pure lattice SU (N) gauge theories with N = S, 10, 12, as the temperature T crosses T c . At these 
large values of N, interactions between singlet objects in the confined phase are suppressed, and a 
Hagedorn picture of string proliferation is most attractive. Because the strings proliferate with an 
effective zero mass, one expect m(T) to vanish at the Hagedorn temperature 7# > T c , which we 
aim to estimate. 

2.1 Methodology and lattice calculation 

Although for N > 3 the first order deconfining transition occurs when m(T c ) > (and therefore 
T c < Th), the Hagedorn transition may still be accessible. This is the case since the deconfining 
transition is strongly first order at larger ,/V, and one can try to use its metastability to carry out 
calculations in the confining phase for T > T c . The way m(T) drops with T tells us how does the 
curvature of the confining minimum in the effective loop potential vanishes. The simplest possi- 
bility is that the point where it vanishes is the spinodal point of the potential, where the confined 
vacuum becomes unstable (and where the barrier between the confining and deconfining vacua dis- 
appears). It is, however, quite possible that 7# does not coincide with the spinodal point, and if the 
latter occurs below 7#, then it may (but need not) interfere with our determination of 7#. This is 
a significant ambiguity that we cannot resolve in the present calculations but the reader should be 
aware of its existence. 

We therefore begin deep in the confined phase and increase T to T > T c , calculating the de- 
crease in m{T), and extrapolating to m = 0. We interpret the result of the extrapolation as the 
Hagedorn temperature, Th- Nevertheless, since we work with finite values of N and volume V, 
tunneling probably occurs somewhere below Th- These tunneling effects and the fact that as m(T) 
decreases, finite volume effects become important, can make an apriori fit for the critical exponents 
unreliable. As a result we first perform fits where we fix the functional behaviour of the temporal 
loop mass to be m(T) = A ■ (Th/T c — T /T c ) v , where V = 0.6715(3) corresponding to three dimen- 
sional XY model (3DXY) or v = 0.5 corresponding to mean field. In addition we also perform fits 
where the exponent v is a free parameter, constrained to be positive and smaller than one. 

Our lattice calculations practicalities are the same as those in Section |TJ We work with L( — 5, 
and L s = 12. In addition, every five sweeps, we measure the correlations functions of improved 
operators for Polyakov lines. The masses are extracted with a variational technique, and the re- 
sults are given in [JO]], where the full detailed description of our lattice study (including statistics, 
thermalization details, choices of initial configuration, monte-carlo results, etc.) is given as well. 
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The physical scale was fixed using the interpolation for the string tension given in [Q] in the 
case of SU (8). For N = 10, 12 we extrapolate the parameters of the scaling function of [Q] in 1/N 2 
from their values at N = 6, 8. In addition we measured the string tension for N = 10 at j3 = 68.80 
on an 8 4 lattice, and for N = 12 at j8 = 99.2, 100.0, on an 8 4 , and a 10 4 lattice respectively. We find 
that the measured values deviate at most by 1.6 sigma from the calculated ones. 

Finally, to check the effect of finite volume corrections, we perform several additional calcu- 
lations of the correlation lengths for SU (10) with L s = 14, and for SU (12) with L s = 16. We find 
that the finite volume effects of the extracted masses are quite small, and are at most on the level of 



1.9 sigma (and mostly lower than that, see [11]). Comparing with the situation close to the second 
order phase transition of the St/ (2) group we find that for N = 2, finite volume effects are 
much more important than for ,/V = 12, which is consistent with standard theoretical arguments that 
predict smaller volume corrections for gauge theories with larger values of N. 

We estimated T /T c by , where (ay/a) c is the value of a\fd at j8 = j8 c on an L t = 5 

lattice 2 . The latter we extract from T c / \[c, which is 0.5819(41) for N = 8. The corresponding value 
for the N = 10, 12 is found by extrapolating in 1 /N according to measured values of T C J \fd(L t = 
5) for /V = 4, 6, 8. This gives 0.5758 for SU (10) and 0.5735 for 51/(12), with an error of about 1%. 



2.2 Results 

Here we focus only on the results of the SU (12) study. We choose to plot the effective string 
tension o e a(T) = m(T) T, in units of the zero temperature string tension Ob, in Fig. |2[ 
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Figure 2: Effective string tension for SU (12) obtained in units of the zero temperature string tension. 
We find that the fits have a rather good £ 2 /dof = 0.30,0.50 when fitting with v = 0.6715,0.5, 
and give Tjj/T c = 1.116(9) — 1.092(6), respectively. The result of the fit with v as a free parameter 
is Vf re e = 0.69, and T H /T C = 1.119, with % 2 /dof = 0.37 (the error here is correlated with the error 



of V and A. (For the iso-surfaces that correspond to confidence levels see [|11|]). With this weak 
preference of v = 0.6715 we cannot unambiguously determine the exponent v (the error on Vf ree is 
quite large, which can be seen from the goodness of all these three fits). Nonetheless, in the case 



of SU (10) we find a stronger preference for v = 0.6715, see [11]. 

The limited statistics prevent us from making statements about the function Tjj(N). This is 
unfortunate, since it is of interest to know how far is Tr/T c from 1 at N = °°. Nevertheless for 



2 This assumes that T c /^fo varies at most very weakly with a(/3), which is in fact what one observes 
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N = 8, 10, 12 we obtain fitted values of T H /y/o ~ 0.62 - 0.68, which is lower than T c j^fo ~ 0.7 
of SU (2) where the phase transition is second order, and therefore may be Hagedorn, T c = Th, or 
provides a lower bound on 7#. Nonetheless, when we plot the effective string tension for N = 10, 12 
as a function of T/^G we find that they are quite similar, but differ from the plot of G e ff(T) of 
N = 2, which may question the identification of deconfinement at N = 2 with a Hagedorn transition. 
Finally, an additional outcome of this work is to confirm that at T = T c the mass of the timelike flux 
loop that couples to Polyakov loops is far from zero at large-Af ((T e ff/ob — 0.5 for N = 8, 10, 12, as 
seen in Fig. ||here and in [|ll|]), which confirms that the transition is indeed strongly first-order. 
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